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ABSTRACT 

We present spatially and spectrally resolved Atacama Large Millimeter/submillimeter Array 
(ALMA) observations of gas and dust in the disk orbiting the pre-main sequence binary AK Sco. 
By forward-modeling the disk velocity field traced by CO J=2—1 line emission, we infer the mass 
of the central binary, = 2.49 ± 0.10 Mq, a new dynamical measurement that is independent of 
stellar evolutionary models. Assuming the disk and binary are co-planar within ^2°, this disk-based 
binary mass measurement is in excellent agreement with constraints from radial velocity monitoring of 
the combined stellar spectra. These ALMA results are also compared with the standard approach of 
estimating masses from the location of the binary in the Hertzsprung-Russell diagram, using several 
common pre-main sequence model grids. These models predict stellar masses that are marginally 
consistent with our dynamical measurement (at ^2a), but are systematically high (by ^10%). These 
same models consistently predict an age of 18 ± 1 Myr for AK Sco, in line with its membership in 
the Upper Centaurus-Lupus association but surprisingly old for it to still host a gas-rich disk. As 
ALMA accumulates comparable data for large samples of pre-main sequence stars, the methodology 
employed here to extract a dynamical mass from the disk rotation curve should prove extraordinarily 
useful for efforts to characterize the fundamental parameters of early stellar evolution. 

Subject headings: protoplanetary disks - stars: fundamental parameters - stars: pre-main sequence 
- stars: individual (AK Sco) 

1. INTRODUCTION 

Precise measurements of the physical properties of pre- 
main sequence (pre-MS) stars are fundamental to testing 
the theoretical predictions of stellar evolution models. 

Such models are in turn the basis for deriving a variety 
of interesting quantities in star-formation research, in¬ 
cluding the masses and ages of individual pre-MS stars 
from secondary properties (such as effective temperature 
and luminosity), the initial mass function (IMF) of star¬ 
forming regions, and the timescale for circumstellar disk 
evolution and planet formation, among many others. To 
successfully test and help refine theoretical models of pre- 
MS evolution, we require a sample of “benchmark” sys¬ 
tems where the fundamental properties are known. The 
key constraint for such systems is a direct (i.e., dynami¬ 
cal) measurement of the stellar mass. 

Unfortunately, few such benchmarks exist, particularly 
at the low end of the mass spectrum {M^ < 2 Mq) where 
the models are in greatest need of calibration (see, e.g., 

Hillenbrand & White 2004; Mathieu et al. 2007; Stas- 
sun et al. 2004, 2006, 2007, 2008; Gennaro et al. 2012). 

Eclipsing binaries (EBs) have long served as important 
empirical touchstones for testing stellar models. How¬ 
ever, the recent review of pre-MS benchmarks by Stas- 
sun et al. (2014) identified only 21 low-mass EBs that 
have sufficiently precise measurements of their physical 
parameters to be suitable for testing evolutionary mod¬ 
els. The same review also presented new evidence that 
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many of the benchmark pre-MS EBs may have their tem¬ 
peratures and/or radii altered by the influence of tertiary 
companions which, while representing interesting physics 
in their own right, render them less suitable to direct 
tests of the evolutionary models that do not include such 
effects (see also Gomez Maqueo Ghew et al. 2012). 

Therefore, additional pre-MS stellar mass benchmarks 
are crucial. Ideally, these would include single stars, or 
binaries without the potentially complicating effects of 
a tertiary component. In particular, single or binary 
pre-MS stars with circumstellar or circumbinary disks 
offer the opportunity to dynamically measure the stellar 
masses using the Keplerian orbital motion (i.e., rotation 
curve) of the gas disk (e.g. Simon et al. 2000; Schaefer 
et al. 2009; Rosenfeld et al. 2012b). 

To be sure, there are additional complications with the 
analysis of such systems, including accretion effects, iso¬ 
lating the individual component masses in the case of a 
close binary (since the dynamical mass is then the to¬ 
tal binary mass), and a systematic (linear) dependence 
on the (still uncertain) distance. However, such sys¬ 
tems are much more common than the intrinsically rare 
EBs, and therefore represent an important opportunity 
to significantly expand the sample of benchmark pre-MS 
stars. This is particularly compelling now, as ALMA has 
started to provide unprecedented sensitivity to molecu¬ 
lar line emission from disks, and the upcoming GAIA 
mission is poised to give accurate distances to these sys¬ 
tems. Just as important as the ongoing determination 
of empirical masses for more of these pre-MS systems 
is the development of robust procedures for their careful 
analysis, including sophisticated molecular line modeling 
approaches and the incorporation of state-of-the-art sta¬ 
tistical methods for comparing the stellar properties to 
current theoretical pre-MS evolution models. 
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Figure 1. (a) Synthesized image of the 1.3 mm dust continuum emission. Contours start at 3 cr and increase by factors of two. (b) The 
0th moment map (velocity-integrated emission; contours) overlaid on the 1st moment map (intensity-weighted velocities; color scale) for 
the CO J=2—l emission. Contour levels start at 3 cr and increase in 10 cr intervals, (c) and (d) The same as (b) for the ^^CO and C^®0 
J=2—1 emission, respectively. Contours are spaced at 3 cr intervals. The synthesized beam is shown in the lower left of each panel. 


AK Sco represents a good case study for this purpose. 
AK Sco is a bright {V « 8.9), pre-MS double-lined spec¬ 
troscopic binary (Andersen et al. 1989) associated with 
the nearby Upper Centaurus-Lupus star-forming region 
(Pecaut et al. 2012, d ^ 140 pc). It is actively accret¬ 
ing and has a massive circumbinary dust disk (Alencar 
et al. 2003). An orbital solution from long-term radial 
velocity monitoring of optical spectra identifies a short 
period (^13.6 days), eccentric (e = 0.47), nearly equal- 
mass pair of F5 stars (Andersen et al. 1989; Alencar et 
al. 2003). Alencar et al. (2003) presented the most com¬ 
prehensive analysis of the system to date, including an 
attempt to determine the stellar masses by constraining 
the (unknown) orbital inclination using geometric argu¬ 
ments (the stars are not known to eclipse) together with 
a simple model of the spectral energy distribution (SED). 
Recently, Anthonioz et al. (2015) resolved the inner edge 
of the AK Sco circumbinary disk with the VLT interfer¬ 
ometer. However, until now there has not been a dynam¬ 
ical measurement of the disk rotation curve with which 
to directly measure the binary mass. 

Here we report the direct determination of the AK Sco 
binary mass through the dynamical measurement of its 
circumbinary disk rotation profile using new data from 
ALMA. Section 2 presents the ALMA observations and 
their calibration. Section 3 covers the methodology for 
modeling the disk gas rotation profile, and thereby the 
binary mass. Section 4 discusses these results together 
with an analysis of other basic stellar properties in the 
context of pre-MS evolutionary model predictions. 

2. OBSERVATIONS AND DATA REDUCTION 

AK Sco was observed using 32 ALMA antennas and 
the Band 6 receivers on 2014 April 10 as part of pro¬ 
gram 2012.1.00496.S (PI Andrews), with baseline lengths 
of 15-350 m. The ALMA correlator was configured to 
process four dual-polarization spectral windows. Two 
windows, covering the CO (230.538 GHz) and 
(220.399 GHz) J=2—1 transitions, spanned a 234 MHz 
bandwidth with 61 kHz (80 m s“^) channels. A third win¬ 
dow sampled G^^O J=2—1 (219.560 GHz) at twice that 
channel spacing (122kHz, or 160m s“^). The remaining 
window used a coarse frequency resolution (128 channels 


of 15.628 MHz width) to probe the 232 GHz (1.3 mm) 
continuum. The observations alternated between AK 
Sco and the quasar J1709—3525 (3° separation) on a ^7 
minute cycle. Some additional brief observations of the 
bright quasar J1626—2951 and Titan were made for cal¬ 
ibration purposes. The observing block lasted ^1 hour, 
with half the time devoted to AK Sco. 

The raw visibility data were calibrated using the fa¬ 
cility software package CASA (v4.2). After applying the 
standard system temperature and water vapor radiome¬ 
ter corrections, the intrinsic passband shape was de¬ 
termined using the observations of J1626—2951 and re¬ 
moved. Gomplex gain variations due to instrumental and 
atmospheric effects were calibrated based on the regu¬ 
lar monitoring of J1709—3525, and the overall amplitude 
scale was set using the observations of Titan. After a sin¬ 
gle iteration of phase-only self-calibration, the reduced 
visibilities were time-averaged into 30 s intervals. The 
local continuum level was subtracted from the spectral 
windows containing emission lines of interest. 

Gontinuum and spectral line emission from all three 
of the targeted transitions were firmly detected. The 
calibrated visibilities were Fourier inverted with natural 
weighting, deconvolved with the CLEAN algorithm, and 
subsequently restored with a FWHM = 0775 x 0760 syn¬ 
thesized beam (with P.A. = 100°). Some representative 
data products are shown together in Figure 1. 

The synthesized map of the 1.3 mm continuum emis¬ 
sion (Fig. la) has an RMS noise level of 45 /iJy beam“^, 
and shows a marginally resolved morphology: an ellipti¬ 
cal Gaussian fit to the visibilities indicates a total flux 
density of 32.65 ± 0.07 mJy (a peak S/N « 400; an ad¬ 
ditional systematic uncertainty of ^10% is imposed by 
the absolute accuracy of the Titan emission model), with 
FWHM dimensions of 0738(±0.01) x 07l2(±0.01) ori¬ 
ented at a position angle of 49 ± 1°. The GO J=2—1 line 
is detected over a velocity span of 24 km s ^ (about 300 
channels at the native resolution). It spans a diameter 
of ^2", and has an integrated intensity of2.21±0.01Jy 
km s“^. For computational simplicity, we average these 
data into 60 channels of 305 kHz (0.4km s“^) width for 
further analysis. The RMS noise level at that resolution 
is 4 mJy beam“^; the peak brightness is 170 mJy beam“^ 
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(peak S/N « 45 per beam and channel). With the same 
averaging, the and isotopologue transitions 

are also detected (albeit over a smaller velocity range 
and at lower significance). Their integrated intensities 
are 0.54±0.01 and 0.20±0.01 Jy km s“^ (with peak S/N 
~ 12 and 5 per beam and channel), respectively. 

3. CO MODELING AND RESULTS 

Our primary goal is to use the observed spectral visi¬ 
bilities that trace the CO J=2—1 emission line to quanti¬ 
tatively characterize the velocity field of the AK Sco gas 
disk, and thereby to measure the total mass of the central 
binary host. We employed a forward-mo deling approach 
that uses a parametric prescription for the disk structure 
(densities, temperatures, and dynamics), and then sim¬ 
ulates the observed visibilities by assuming the molecu¬ 
lar level populations are in local thermodynamic equilib¬ 
rium (LTE) and propagating synthetic photons through 
the model structure with the RADMC-3D radiative transfer 
code.^ Previous work with similar intentions (Simon et 
al. 2000; Pietu et al. 2007; Schaefer et al. 2009; Guilloteau 
et al. 2014) has shown that resolved measurements of the 
disk rotation curve can constrain the central stellar mass 
(M^) with high precision (^few percent, although with a 
linear dependence on distance). However, in practice this 
approach is complicated: it involves radiative transfer 
modeling with a large number of (unrelated) disk struc¬ 
ture parameters. Given that complexity, it is important 
to develop tests to validate its absolute accuracy. Rosen- 
feld et al. (2012b) used the circumbinary disk around 
the V4046 Sgr system to demonstrate that this approach 
is robust, in that it provides an estimate consistent 
with the constraints from spectroscopic monitoring of the 
stellar radial velocities. The AK Sco system provides an¬ 
other rare opportunity to test the methodology. 

We adopted the two-dimensional (axisymmetric) para¬ 
metric disk structure model that is described in detail by 
Rosenfeld et al. (2012b). The temperature structure is 
vertically isothermal, and has a power-law radial distri¬ 
bution with index q and a normalization at 10 AU (Tio). 
The radial surface density profile is the standard Lynden- 
Bell & Pringle (1974) similarity solution for a viscous ac¬ 
cretion disk, essentially a power-law with an exponential 
taper at large radii. The gradient parameter 7 is fixed 
to unity, so the profile is described by a characteristic ra¬ 
dius Vc and a total GO mass Mco- The densities are dis¬ 
tributed vertically under the assumption of hydrostatic 
equilibrium. The bulk velocity field of the gas is assumed 
to be Keplerian, and dominated by the binary mass (as¬ 
suming Mdisk/AT* ^ 1). The line width is calculated 
as the quadrature sum of thermal and non-thermal (i.e., 
turbulent) broadening terms, with the latter denoted as 
a constant velocity width The physical structure of 
the model is fully characterized by six free parameters, 

^disk — {^lO? Q.-) ^co-) 'f'ci 

For a given set of these structure parameters, the 
RADMC-3D radiative transfer code was used to cal¬ 
culate the corresponding molecular excitation levels 
(appropriately assuming LTE for this transition; see 
Pavlyuchenkov et al. 2007) and ray-trace synthetic model 
spectral images at high resolution (2.5 AU pixels). To do 
so, we specified an additional set of parameters. The 

^ ita.uni-heidelberg.de/~dullemond/software/radmc-3d/ 


disk inclination id is defined with respect to the rotation 
axis, such that id = 0° corresponds to a face-on viewing 
geometry with the rotation axis pointing toward the ob¬ 
server (and the disk rotating counterclockwise), id = 90° 
is an edge-on orientation, and id = 180° is again face- 
on but with the rotation axis pointing away from the 
observer (and therefore an apparent clockwise sense of 
rotation). The position angle ip represents the projec¬ 
tion of the angular momentum vector of the disk onto 
the sky (as typical, oriented E of N). The disk center 
is characterized by offsets (Jq,, 5s) relative to the phase 
center. Along with the distance d and systemic velocity 
i;sys, there are six additional free parameters, ^obs = 

^6-) '^sys: 

The end result is a set of high resolution model im¬ 
ages that specify the sky-projected GO J=2—1 surface 
brightness as a function of position and frequency for 
any set of parameters 6 = {^disk, ^obs}- We then em¬ 
ployed the FFTW algorithm (Frigo & Johnson 2005) to 
Fourier transform these spectral images, and then per¬ 
formed a band-limited interpolation of the complex visi¬ 
bilities onto the same spatial frequencies (ix, v) sampled 
with ALMA (using the spheroidal gridding functions ad¬ 
vocated by Schwab 1984) to acquire a set of model visi¬ 
bilities Mu^y{0). To quantify the model quality given the 
data Du^y^ we computed a simple likelihood function 

p{D\d) = '[lw^,,\D^,,-MuAOp, ( 1 ) 

u,v 

where Wu,v are the visibility weights (determined from 
the radiometer equation). Eq. (1) is identical to a log- 
likelihood function described by the sum of the standard 
values (real and imaginary) over all the observed spa¬ 
tial frequencies and velocity channels. To explore the 
posterior distribution of the model parameters, we em¬ 
ployed the Markov Ghain Monte Garlo (MGMG) tech¬ 
nique with a Metropolis-Hastings (M-H) jump proposal. 

We assumed uniform priors in all parameters except 
d, since the ALMA data do not constrain the distance. 
A trigonometric parallax for AK Sco is available from 
Hipparcos^ although the considerable optical variability 
is problematic: the original catalog has d = I44I25 PC 
(Perryman et al. 1997), but the van Leeuwen (2007) re¬ 
vision argues for a closer value (d = 102^^7pc). Pecaut 
et al. (2012) measured a kinematic parallax correspond¬ 
ing to d = 144 ± 12 pc. Recently, Anthonioz et al. (2015) 
combined astrometric and radial velocity monitoring to 
estimate d = 141 ± 7 pc. We adopted a Gaussian prior 
on d based on the weighted average of these latter two 
measurements, with mean 142 pc and cr = 6pc. 

After a few preliminary chains were run, we empiri¬ 
cally tuned the covariance of the M-H jumps to match 
the morphology of the posterior distribution, yielding a 
more efficient exploration of parameter space. To de¬ 
termine the final posterior distributions, we ran multiple 
independent chains using different initializations. After a 
conservative period of burn-in, we computed the Gelman- 
Rubin statistic (Gelman et al. 2013) for each parameter 
over the ensemble of chains to assess convergence.^ 

® The code used to perform the analysis described here is 
open source and freely available under an MIT license at https: 
//github.com/iancze/JudithExcalibur. 
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Figure 2. Channel maps at 305kHz (0.4km s resolution showing the CO J=2—l data, the corresponding model, and the residuals 
(the latter two imaged in the same way as the data). Contours are drawn at 3cr (12mJy beam“^) intervals, and the synthesized beam 
dimensions are shown in the lower left corners of each set of channel maps. The LSRK velocities are labeled in the top left of each panel. 
The color within each channel map corresponds to the velocity sampled in the moment maps (Figure 1). 
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Table 1 

Inferred Parameters for AK Sco 


Parameter 

Value 

Parameter 

Value 

Tio (K) 

92 ± 4 

(°) 

109.4 ± 0.5 

q 

0.51 ± O.OI 

¥^(°) 

I4I.I ± 0.3 

logMco (Me) 

-0.65 ± 0.22 

Vsys^ (km s 

+5.49 ± 0.06 

rc (AU) 

14.3 ± 1.2 

5a (") 

0.053 ± 0.002 

-M* (-Mq) 

2.49 ± O.IO 

5« (") 

0.045 ± 0.002 

^ (km s-1) 

0.31 ± 0.02 

d (pc) 

143.6 ± 5.7 


Note. — The quoted values represent the “best-fit”, the peaks 
of the marginal posterior distributions. The uncertainties corre¬ 
spond to the 68.3% (~1 cr) confidence intervals. 

^ In the LSRK frame, for the standard radio definition. The 
corresponding heliocentric value is —1.92 ± 0.06, consistent with 
the Alencar et al. (2003) value derived from optical spectroscopy. 

The parameter values inferred from modeling the ob¬ 
served spectral visibilities binned to 305 kHz (0.4 km s“^) 
resolution are listed in Table 1. A trial fit of the data 
at the best available (quasi-independent) spectral reso¬ 
lution of 122 kHz (0.16 km s“^) recovered these same val¬ 
ues. The disk structure parameters are typical for similar 
Class H systems, although the characteristic radius and 
mass are on the small side. Assuming a standard CO/H2 
abundance ratio (^10“^), the total gas mass would be 
^0.007 Mq - roughly consistent with what would be in¬ 
ferred from the continuum emission (e.g., Andrews et 
al. 2013). The inferred total mass of the central binary 
is 2.49 ± 0.10 Mq, where the uncertainty is dominated 
by the distance prior. At a fixed d, the constraint can 
be framed as M^/d = 0.01731 ± 0.00015 Mq pc“^; i.e., 
the formal precision on the stellar mass is ^1% for a 6- 
function prior on d. Figure 2 shows a direct comparison 
of the data and best-fit model in the image plane, demon¬ 
strating (through the absence of significant residuals) the 
fit quality. The reduced of the best-fit model is 1.08. 

4. DISCUSSION 

We have presented ALMA Cycle 1 observations around 
a wavelength of 1.3 mm that resolve the dust and molec¬ 
ular line emission from the disk that orbits the double- 
lined spectroscopic binary AK Sco. To our knowledge, 
these are the first observations that confirm the pres¬ 
ence of a substantial molecular gas reservoir in this sys¬ 
tem. The disk itself is relatively small (a characteris¬ 
tic radius of ^14 AU), but still contains a modest mass 
(^5-10 Mjup in total, for standard assumptions about 
opacities and abundances). We have focused on a de¬ 
tailed modeling analysis of the spatially and spectrally 
resolved CO J=2—1 line emission, to map out the disk 
velocity field and make a dynamical estimate of the to¬ 
tal mass of the host binary. Adopting a well-motivated 
prior on the distance (142 ±6 pc; see Sect. 3), we inferred 
a precise (4%) combined stellar mass estimate given the 
ALMA spectral visibilities, = 2.49 ± 0.10 Mq. 

4.1. Comparison to Previous Results 

There are two independent dynamical constraints on 
the AK Sco binary mass available in the literature. The 
first comes from an optical spectroscopic campaign to 
monitor the radial velocity (RV) variations of the double- 
lined system, performed by Alencar et al. (2003). They 
measured a precise orbital period (^13.6 d) and mass 



Figure 3. The (marginal) joint posterior probability distributions 
for {M*, inclination} and For reference, we have overlaid 

the nominal parameter degeneracies (±I cr) that accompany each 
measurement in isolation, assuming no external constraints about 
either parameter but precise knowledge of all other parameters. 

ratio (0.987 ± 0.005), a substantial eccentricity (e^ = 
0.471 ± 0.002) and a total mass constraint sin^ A = 
2.114 ± 0.010 Mq where A is the (unknown) inclination 
of the binary orbit. These measurements do not de¬ 
pend on d. If we assume the gas disk and binary or¬ 
bits are co-planar (A = id)^ then our constraint on id 
would convert the RV constraint on the binary mass 
to 2.52 ± 0.03 Mq, identical (within 1%) to our disk- 
based dynamical mass estimate. These two dynamical 
mass constraints remain consistent (at the 1 a level) so 
long as the binary and disk orbital planes are aligned 
within ^2°, compatible with the statistical constraints 
on the mutual inclination distributions found for exo¬ 
planet systems (Tremaine & Dong 2012; Figueira et al. 
2012; Fabrycky et al. 2014). Additionally, by combin¬ 
ing the distance-independent measurements of /d and 

sin^ A, we derive a dynamical distance to the system 
of d = 145.5 ± 2.0 pc, which is consistent with our initial 
assumptions in the form of the distance prior. 

This consistency between mass constraints, along with 
the similar results for the V4046 Sgr binary (Rosenfeld et 
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al. 2012b), effectively validates the absolute accuracy of 
the disk-based dynamical mass measurement technique. 
Whereas the V4046 Sgr measurement involved a large 
circumstellar disk {vc = 45 AU), we have demonstrated 
that—given the sensitivity and resolution of ALMA—the 
dynamical mass technique also works for much smaller 
disks like AK Sco {vc = 14 AU). In the coming age of 
ALMA, this should be the workhorse approach for de¬ 
termining masses for large, statistically relevant, samples 
of young stars. Most importantly, it is the only method 
capable of performing that task for single stars. 

The second mass constraint comes from the recent 
work by Anthonioz et al. (2015), which combined multi¬ 
epoch iL-band interferometric observations with the 
Alencar et al. (2003) RV data. In the context of a sim¬ 
ple model that includes the binary and a narrow dust 
ring, these astrometric constraints suggest inclinations 
ofi^ = 115±3° and iring = 121 ±8°, respectively, which 
are ^1.8 and 1.4 <j larger than the id we infer from the 
ALMA data. In conjunction with the RV constraints, 
this larger orbital inclination suggests a higher binary 
mass, = 2.80 ± 0.11 Mq, a 2 a discrepancy with re¬ 
spect to the disk-based mass estimate. For reference. 
Figure 3 shows the mass constraints as a function of both 
inclination and distance. 

There are various scenarios that might alleviate this 
(modest) tension. The difference between i^ing and 
id might be explained with a warp in the disk struc¬ 
ture, which would perturb the line-of-sight projection of 
the velocity field measured in the CO line and thereby 
accommodate a higher estimate from the ALMA 
data (e.g., Rosenfeld et al. 2012a; Rosenfeld et al. 2014; 
Marino et al. 2015). Given the eccentricity of the bi¬ 
nary, one might also expect that the inclination esti¬ 
mates could be biased due to the assumptions of circular 
disk and ring structures. For the Anthonioz et al. re¬ 
sults, an imperfect ring model can affect the estimate of 
i* since the binary itself is not well-resolved (the semi¬ 
major axis of the apparent orbit is only 1 mas, about one 
third of the resolution). We also explored an eccentric 
disk model, where the structure is built up from a series 
of infinitesimal apse-aligned elliptical rings: the ALMA 
data rule out a significant (mean) eccentricity in the gas 
disk, with Cd < 0.04 (at 99.7% confidence). In the end, it 
may just be that the (admittedly) simplistic ring model 
adopted in the preliminary analysis of Anthonioz et al. 
(2015) can be improved (the quoted reduced of their 
fit is ^2), and such modifications would bring the results 
into agreement. Ultimately, it would be interesting to 
combine all the measurements in a joint analysis. 

4.2. Comparison to Pre-MS Evolution Models 

To compare the disk-based dynamical mass with pre¬ 
dictions from pre-MS evolution models, we performed the 
standard analysis of estimating stellar parameters from 
the Hertzsprung-Russell (H-R) diagram. The UBVRI 
photometry compilation of Jensen & Mathieu (1997) was 
used to construct the AK Sco SED. The near-infrared is 
contaminated by the dust disk, and so was excluded from 
our analysis. A base parametric SED model was con¬ 
structed from synthetic photometry in the BT-Settl 
catalog of stellar models (Allard et al. 2003), interpo¬ 
lated for any given effective temperature (Teff) and sur¬ 
face gravity (log^) assuming a fixed solar metallicity. 


Table 2 

Evolutionary Model Predictions 


Model Grid 


[Me] 

r* 

[Myr] 

Siess et al (2000) 

2.70 

± 

0.06 

19 


3 

Tognelli et al. (2011) 

2.71 

± 

0.07 

18 


2 

Dotter et al. (2008) 

2.60 

± 

0.07 

17 


2 

Baraffe et al. (2015) 

2.68 

± 

0.05 

19 


2 


Note. — The quoted uncertainties corre¬ 
spond to the 68.3% (~1 a) confidence intervals. 


This base model was adjusted for extinction using the 
Eitzpatrick (1999) reddening curve with Ry = 4.3 (as 
determined by Manset et al. (2005)) and then scaled by 
the squared ratio of the radius (R*) and distance (d). 
We assumed the components of the AK Sco binary are 
identical. 

This SED model has five free parameters, ^sed = {Feff, 
log^. Ay, R^^ d}, but the photometric data alone cannot 
uniquely constrain all of them. To aid in the inference of 
the stellar properties, we imposed simple Gaussian pri¬ 
ors on d (as in Sect. 3) and Teff. The latter was based 
on the E5 spectral classification and de-reddened color 
indices (Andersen et al. 1989), which we associate with 
Teff = 6450T150K (e.g., Bessell 1979; Popper 1980; Gray 
1992; Gasagrande et al. 2010; Pecaut & Mamajek 2013). 
AK Sco is highly variable, with erratic changes much 
larger than the uncertainties on individual measurements 
(e.g., Andersen et al. 1989; Alencar et al. 2003). To deal 
with that variability, we incorporated a nuisance “jitter” 
parameter (a) at each band that characterizes the ad¬ 
ditional dispersion (assuming a Gaussian distribution). 
With this parametric model setup, we explored the pos¬ 
terior distribution of 6 = {^sed, crp, crp, cry, ct/} 
conditioned on the SED data using MGMG with the en¬ 
semble sampler emcee (Eoreman-Mackey et al. 2013). 

The modeling results are shown in Eigure 4. We found 
Teff = 6365T155K, log^^ = 3.5T0.5, R* = 1.43±0.07Ro, 
and Ay = 0.70 ± 0.1 mag for each star (and d = 142 ± 
6pc, as expected given the prior). The corresponding 
luminosity of each component is = 3.0 ± 0.5 Tq. The 
variability dispersion terms, in terms of a flux density 
fraction, range from 0.1 (/-band) to 0.2 (U-band). 

Models of pre-MS evolution predict the joint behav¬ 
ior of Teff and R^ (or equivalently T*) as a function of 
and age (r*). Following the Bayesian formalism of 
Jprgensen & Lindegren (2005), we mapped the posterior 
constraints from the SED modeling into an inference on 
{M^, r*} for several evolutionary model grids (for some 
practical examples, see Rosenfeld et al. 2012b; Andrews 
et al. 2013). This approach is also illustrated in Fig¬ 
ure 4. The AK Sco binary mass and age inferred from 
these model grids (Siess et al. 2000; Dotter et al. 2008; 
Tognelli et al. 2011; Baraffe et al. 2015) are listed in Ta¬ 
ble 2. Their weighted means and standard deviations are 
M^ = 2.68 ± 0.03 Mq and = 18±lMyr. These model- 
predicted masses are consistent, but are all ^1.5-2 a 
higher than the dynamical mass inferred from the ALMA 
data. The ages are also consistent with the age of the 
Upper Centaurus-Lupus association ((r*) = 16 ± 1 Myr; 
Pecaut et al. 2012). 

The 7-10% discrepancy between the pre-MS model 
masses and the disk-based dynamical mass is slightly 
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Figure 4. top: The SED compared with the best-fit model. The 
residuals are shown overlaid with the inferred “jitter” terms in 
each band to account for the systematic variance from variability. 
middle: The Teff version of the H-R diagram, with the marginal 
posteriors inferred from the SED modeling shown as 1 and 2 a 
contours. The Dotter et al. (2008) pre-MS model mass tracks are 
overlaid, with color scales indicating different ages (note that these 
values refer to each individual star in the AK Sco binary), bottom: 
The joint mass and age constraints from various pre-MS models, 
shown as 1 (7 contours and marginalized distributions. The gray 
band marks the disk-based constraint on M* (1 cr). The Baraffe et 
al. (2015) models do not explore the parameter space M* > 1.4Mq. 


young EBs in this mass range (Stassun et al. 2014). We 
suspect this modest mismatch might be attributed to the 
complexity of the AK Sco binary environment. The AK 
Sco stars are separated by only r^llR^ at periastron (ev¬ 
ery r^2 weeks), which leads to accretion bursts onto the 
stellar surface (Gomez de Castro et al. 2013). Perhaps 
related, the AK Sco stars have unusually broad ultravi¬ 
olet lines that modulate with the binary period and are 
indicative of perturbed, hot (^60,000 K) “atmospheres” 
that extend out to 5 (Gomez de Castro 2009). This 
behavior is compounded by the erratic variability noted 
in broadband photometry and the Balmer emission lines 
(Alencar et al. 2003), suggesting additional complexity in 
the accretion/outflow behavior. This dynamic, compli¬ 
cated environment likely impacts the physical structures 
and evolution of the stars, and is (obviously) not included 
in the pre-MS models. Alternatively, and perhaps more 
likely, the standard means of estimating parameters like 
Teff in such a situation should also be affected: a small 
shift (< 100 K cooler) would bring the measurements into 
agreement. 

5. SUMMARY AND CONCLUSIONS 

We have analyzed ALMA observations of the CO 
J=2—1 transition from the AK Sco circumbinary disk. 
The main conclusions of this work include: 

• A relatively compact disk in orbit around the AK Sco 
binary, with ^5-10 Mjup of gas and dust, is detected in 
the 1.3 mm continuum and main isotopologues of CO 
(the CO, ^^CO, and C^^O J=2—1 transitions). This 
suggests an unexpectedly long-lived (^ISMyr) disk of 
primordial origin, as opposed to a second-generation 
debris disk. 

• We employed a parametric disk structure model front- 
end fed into the radiative transfer code RADMC-3D to 
generate synthetic spectral visibilities to compare with 
the observations in the MCMC framework. The related 
software is provided as an open source resource to the 
community. The results offer a high-quality dynamical 
measurement of the binary mass, = 2.49 ± 0.10 Mq, 
that is independent of pre-MS evolution models. 

• This disk-based dynamical mass estimate is in good 
agreement with the constraints from radial velocity 
monitoring of the binary, so long as the disk and binary 
orbital planes are aligned within ^2°. There is minor 
tension with a recent combined astrometric + radial 
velocity analysis, although we expect the comparison 
could be improved with a consistent joint analysis. 

• With the standard approach of estimating stellar pa¬ 
rameters from the H-R diagram, we make comparisons 
between pre-MS evolutionary model predictions and the 
dynamical mass estimated here. These models suggest 
a slightly higher (by 7-10%) stellar mass for AK Sco, 
in modest disagreement (at 1.5-2 a) with our results: 
this discrepancy may be attributed to the complicated 
accretion and interaction environment of the binary. 
The model-dependent ages (18 ± 1 Myr) are consistent 
with the proposed AK Sco membership in the Upper 
Centaurus Lupus association. 


larger than the typical level of disagreement noted for 















• The overall consistency between estimates for AK 
Sco validates the absolute accuracy of the disk-based dy¬ 
namical mass technique. This method has great promise 
in the ALMA era, since it is uniquely capable of provid¬ 
ing precise (few %) masses of statistically large samples 
of single pre-MS stars that can be used to test and cali¬ 
brate models of early stellar evolution. 

The authors appreciate the helpful advice of Eric Stem- 
pels and Maxwell Moe; help with an eccentric disk model 
from Meredith Hughes, Sam Factor, and Eugene Chiang; 
and the technical support regarding our many ALMA 
calibration questions from Jennifer Donovan Meyer. IC 
is supported by the NSF Graduate Fellowship and the 
Smithsonian Institution. SA acknowledges the very help¬ 
ful support provided by the NRAO Student Observing 
Support program related to the early development of this 
project. This paper makes use of the following ALMA 
data: ADS/JAO.ALMA#2012.1.00496.S. ALMA is a 
partnership of ESO (representing its member states), 
NSF (USA), and NINS (Japan), together with NRC 
(Canada) and NSC and ASIAA (Taiwan), in coopera¬ 
tion with the Republic of Chile. The Joint ALMA Ob¬ 
servatory is operated by ESO, AUI/NRAO, and NAOJ. 
Figure 3 was generated with triangle.py (Foreman- 
Mackey et al. 2014). This research made extensive use of 
the Julia programming language (Bezanson et al. 2012) 
and Astropy (Astropy Collaboration 2013). 
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